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TITLE OF THE INVENTION 

Tomographic Device and Method Therefor 

FIELD OF THE INVENTION 

[0001] 

5 The present invention relates to a tomographic 

device and a method therefor, which produces highly 
accurate tomographic images of a subject from 
projection data obtained from a radioactive ray source 
and radioactive ray detector constituted movable 
10 relative to the subject in the circumferential 
revolving direction and in the body axis direction. 
CONVENTIONAL ART 
[0002] 

Recently, a multi detector row computer tomography 
15 (herein after will be called as MDCT) has appeared, in 
which a plurality of detector rows are arranged in the 
circumferential revolving axis direction. In 
comparison with a single detector row computer 
tomography (herein after will be called as SDCT) , since 
20 the MDCT is provided with detectors in broad width by 
arranging a plurality of detector element rows along 
the circumferential revolving axis direction, a broad 
imaging area can be covered at one time. Further, with 
the MDCT, when the subject is moved relatively with 
25 further higher speed, the scanning time is shortened, 
thereby, artifacts due to body motion such as breathing 
can be reduced and a resolution in the circumferential 



revolving axis direction can be significantly enhanced. 
Fig.l is diagrams showing fundamental constitutional 
differences between the SDCT and MDCT . The SDCT is 
provided with X-ray detectors 11 in a single row with 
5 respect to a single X-ray source 10 as shown in Fig.l 

(A) , and the MDCT is provided with X-ray detectors 12 
in a plurality of rows (in the drawing, 8 rows) with 
respect to a single X-ray source 10 as shown in Fig.l 

(B) . 

10 [0003] 

In the case of MDCT, since every X-ray detector 
rows locate in different slanting angles with respect 
to the circumferential revolving axis direction, 
parameters for specifying the projection data 

15 increase such as channel, row and slanting angle and 
now the image reconstruction methods therefor are 
complicated and diversified. Under this circumstance, 
a variety of image reconstruction algorithms are 
proposed such as 3-dimensional Radon transformation 

20 method and 3-dimensional back projection method 
{3-dimensional reconstruction methods) as a 
reconstruction algorithm when a further accuracy is 
required, and including weighted back pro j ect ion method 
for helical correction (2-dimensional reconstruction 

25 method) used for MDCT requiring high speed computation, 
which is formulated by improving a weighted back 
projection method for helical correction being used for 



the for the SDCT. 
[0004] 

Among these image reconstruction methods, in the 
weighted back projection method for helical correction, 
5 which is a 2-dimensional image reconstruction method, 
the reconstruction time per one tomographic image is 
short in that from a few seconds to a few tens seconds. 
In an actual device, when a dedicated hardware such as 
DSP board and ASIC is used, the images can be 

10 reconstructed in about 0.2-0.5 seconds per one 
tomographic image. Further, an amount of memory 
required for producing projection data equivalent to 
one row detectors from plural row detectors and for 
performing the 2-dimensional back projection is almost 

15 the same as required in the SDCT and is fully 
satisfactory in view of the cost thereof. Accordingly, 
in MDCTs provided with detector rows such as 2 rows and 
4 rows, the improved 2-dimensional reconstruction 
method is generally employed. 

20 [0005] However, since the weighted back pro j ectibn 

method for helical correction uses an algorithm which 
neglects beam slanting (cone angle) of X-rays in the 
circumferential revolving direction, an image quality 
of MDCTs with detectors of more than 16 rows is 

25 significantly deteriorated due to the influence of the 
cone angle, which reduces the diagnostic accuracy of 
tomographic devices. For this reason, the application 
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of the weighted back projection method for helical 
correction is limited to MDCTs with detectors of about 
2^-8 rows of which influence of cone angle is 
comparatively small . 
5 [0006] 

Recently^ since the number of rows of the detectors 
increases, image reconstruction methods with high 
accuracy for MDCTs having a broad cone angle have been 
generally studied- Among these, although 3-dimensional 

10 Radon transformation method is a precise image 
reconstruction method, an extremely long computing time 
such as from a few tens minutes to a few hours for 
obtaining one slice image is required, which prevents 
practical use thereof. 

15 [0007] 

On the other hand, although 3-dimensional baclc 
projection method is an approximal image reconstruction 
method, however, is an image reconstruction method with 
comparatively high accuracy which talces into account 

20 of the cone angle, of which computation time for one 
slice image is about from a few minutes to a few tens 
minutes, and when a dedicated hardware is used, the 
computation time will be further shortened, therefore, 
the method performs a comparatively high speed 

25 computation and practical. For this reason development 
of realizing an MDCT implementing the 3-dimensional 
baclc projection method is being advanced. 
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[0008] 

One of the problems of the image reconstruction 
method with the highly accurate 3-dimensional back 
projection method correctly taking into account of the 
5 cone angle is to increase significantly the memory 
amount required when performing the image 
reconstruction computation in comparison with the 
2-dimensional back projection method used for the 
conventional SDCTs. Namely, in the back projection 

10 computing unit therein, data (pro jection data) required 
for the back projection are read out from the hard disk, 
store the same in a high speed memory (for example, a 
cache memory) and the back projection processing is 
executed with the data in the high speed memory. In this 

15 instance, when the amount of data to be processed is 
large, a part of the data are stored in a low speed large 
capacity memory (for example, DRAM) , and when data 
necessary for the computation do not exist in the high 
speed memory, the data are successively read out from 

20 the low speed memory, the data in the high speed memory 
are renewed and after the renewing, the processing is 
performed. Since the high speed memory is generally 
expensive, the capacity of the high speed memory is 
mostly small in comparison with that of the inexpensive 

25 low speed memory. 
[0009] 

Now, the amount of memory (amount of data to be 
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processed) required for the back projection processing 
is discussed, in the weighted back projection method 
for helical correction using the 2-dimensional back 
projection method, helical correction projection data 
5 for one detector row are produced from a plurality of 
data through interpolation. Since the back projection 
processing is performed for every view, the necessary 
amount of memory (amount of data to be processed) is 
that for one view. Namely, the amount of memory 

10 necessary for one time (for one view) back projection 
is that for one row x that for number of channels. For 
example, when assuming that the number of channels is 
1000 [ch] , the necessary amount of memory is about 2 
[Kbyte] (= 1000 [ch] x 1 [row] x 2 [byte] ) . On the other 

15 hand, in the S-dimensional Radon transformation method 
and the 3-dimensional back projection method, since the 
detector data from a plurality of rows have to be treated 
as they are, the amount of memory necessary for one time 
(for one view) back projection processing increases in 

20 proportion to the number of detector rows. For example, 
in the case of MDCT having detectors of 128 rows, the 
amount of memory necessary for the weighted back 
projection method for helical correction is 128 times 
in that about 256 [Kbyte] . 

25 [0010] 

As indicated above, since the necessary amount of 
memory (amount of data to be processed) increases and 
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the data can not be stored in the high speed memory within 
the processing unit, a memory swapping is required in 
such a manner that the data to be processed are 
temporarily stored in a low speed memory connected 
5 outside of the processing unit and the processing is 
performed while successively replacing data depending 
on the necessity. In this instance, the processing speed 
depends on the data transmission speed between the high 
speed memory and the low speed memory and a processing 

10 speed more than the data transmission speed can not be 
obtained, which causes a delay of the processing time. 
Further, even when a dedicated hardware is prepared, 
like delay depending on the data transmission speed is 
caused. As will be understood from the above, in order 

15 to obtain a processing speed more than the data 
transmission speed, the capacity of the expensive high 
speed memory has to be increased, however, which 
significantly increases the cost thereof and is not 
desirable . 

20 [0011] 

Another problem of the image reconstruction method 
with the highly accurate 3-dimensional back projection 
method correctly dealing the cone angle is to increase 
the processing time in comparison with the 

25 2-dimensional back projection method used for the 
conventional SDCTs. In the 2-dimensional back 
projection method used for such as SDCTs and MDCTs with 
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4 detector rows, since the back projection is performed 
by using scanning data on an imaginary circular orbit 
of a single detector row which are obtained by helically 
correcting scanning data on a helical orbit by weighting, 
5 the address calculation of the detectors in the row 
direction was not necessary. On the other hand, in the 
3-dimensional back projection method, in order to 
access detector data in a plurality of rows, calculation 
of detector address (addressing) in the channel 

10 direction and in the row direction is required through 
complicating computation according to, for example, the 
following equations. Further, the processing equations 
for the addressing applied to the 3-dimensional back 
projection method according to the present invention 

15 are not limited to the following equations (l)-'{6), but 
a variety of processing equations can be applied 
therefor. 

ti(Ki,yi,*)=Xi • cos<i) + y, ' sin<f) • — (1) 
(z.-zc(xt,yu<^))' SID . . . (2) 

T -[(i^ +arcsini — — — — /] 
25 L(xi,y,.(f>)=D(x,.y,.<^.)+Wi(xi,3^,<^») * * ' (4) 



20 



D(Xi,y,,</))=Vsop--tr • • • (5) 
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[0012] 

Herein, xi^yi^zi shows a coordinate position of 
5 voxsel 1 in an image reconstruction area, 0 shows a 
circumferential revolving position of parallel beams, 
w, t, V are coordinate axes of the detectors wherein 
w represents an axis in advancing direction of the 
parallel beams, t represents an axis in perpendicular 

10 direction (the channel direction of the parallel beams) 
to the advancing direction and v represents an axis of 
the detectors in the circumferential revolving axis 
direction, wi, ti, vi represents a coordinate position 
on the w, t, V axes when parallel beams in 0 phase pass 

15 the coordinate position (xi, yi). SID represents the 
distance between the radioactive ray source and the 
detectors and SOD represents the distance between the 
radioactive ray source and the center of rotation. Zg 
represents the position of the radioactive ray source 

20 in z axis direction and Zso represents the position Zg 
when the circumferential revolving phase of the 
radioactive ray source is zero. 
[0013] 

Fig. 2 is a diagram for explaining the general idea 
25 of image reconstruction method in SDCT. Fig. 3 is a 
diagram for explaining the general idea of image 
reconstruction method in MDCT. In the computation with 



.J 
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the image reconstruction method for the MDCT as shown 
in Fig. 3, the addresses of reconstructed image 30 in 
X and y directions vary linearly in the channel 
direction and non-linearly in the row direction. On the 
5 other hand^ in the computation with the image 
reconstruction method for the SDCT as shown in Fig. 2^ 
in response to the linear address change . of 
reconstructed image 20 in x and y directions, the 
corresponding detector addresses change linearly in the 

10 channel direction. As will be understood from the above, 
the 3-dimensional back projection method shows a 
drawback of accompanying a significant delay for 
processing data required for the back projection 
computation due to 2-dimensional addressing and complex 

15 addressing. In particular, the addressing in the row 
direction in the 3-dimensional back projection method 
requires to use extremely complicated non-linear 
functions and simplification by modifying the functions 
are difficult, which is the great cause of the data 

20 processing delay. 
[0014] 

A tomographic device which takes into an account 
of resolving such problems is, for example, disclosed 
in JP-A-2003-24326, in which back projection 
25 computation (tomographic image reconstruction 
calculation) for performing back projection on 
2-dimensional or 3-dimensional tomographic image 
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reconstruction area imaginarily set on a region of 
interest of a subject is performed by a computer for 
every divided region formed by dividing the tomographic 
image reconstruction area, thereby, the tomographic 
5 image reconstruction calculation can be performed 
orderly for every optimum region (the divided region) 
determined in view of the cache memory size, as a result, 
data reuse rate in the cache memory is increased, data 
access times with the memory is reduced, total data 
10 transference time . for the tomographic image 
reconstruction is shortened and the tomographic image 
reconstruction calculation time is shortened. 
[0015] 

However, as disclosed in JP-A-2003^24326, in order 
15 to perform back projection on the divided tomographic 
image reconstruction region stored in the high speed 
cache memory, all of the projection data thereof are 
stored in the cache memory, which increases the amount 
of high speed memory and is not preferable in view of 
20 the cost of the device. Further, in the case of MDCTs 
of which increase in number of detector rows is extreme, 
the amount of data to be processed becomes significant, 
which causes significant increase of the high speed 
memory amount and prevents processing cost reduction 
25 and high speed processing. 

An object of the present invention is to provide 
a tomographic device and a method therefor which 
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suppress capacity increase of a high speed memory, 
prevent significant increase of processing cost and 
permit to produce high quality tomographic images in 
high speed. 
5 SUMMARY OF THE INVENTION 

[0016] 

One feature of a tomographic device according to 
the present invention, which detects penetration light 
penetrated through a subject with a detection means 

10 arranged in 2-dimension and produces 3-dimensional 
tomographic image of a region of interest of the subject 
from the detected projection data, is to provide 
processing means which divides an image reconstruction 
area of the subject into a plurality of image data 

15 segment regions, extracts among the projection data 
detected by the detection means a projection data 
segment region necessary for back projecting on the 
image data segment region and performs 3-dimensional 
back projection computing processing for every image 

20 data segment region by making use of the extracted 
projection data segment region. 

The tomographic device according to the present 
invention is constituted in such a manner that in order 
to reduce a high speed memory amount necessary for the 

25 back projection processing the image reconstruction 
area is divided into a plurality of small regions (image 
data segment region) , a minimum projection data segment 
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region necessary for the back projection computation 
processing is extracted for every divided image data 
segment region among the projection data obtained from 
the imaging and the back projection computing 
p processing is performed for every image data segment 
region of the small region by making use of the data 
of the extracted projection data segment region. 

Further^ in order to maximize the utilization of 
the capacity of the high speed memory, it is preferable 

10 that the size of the projection data in the view 
direction is determined depending on the capacity of 
the high speed memory which can be utilized during the 
image reconstruction processing. Further, in order to 
reduce the complexity of the image reconstruction 

15 processing, it is preferable, when dividing the image 
reconstruction area into image data segments, to divide 
the same into small segments having the same size. Still 
further, in order to reduce the amount of memory 
necessary for data processing at one time, it is 

20 preferable that the size of the projection data segment 
region in the view direction is to be determined one 
for a single view. 
[0017] 

Another feature of the tomographic device 
25 according to the present invention is that in the 
tomographic device provided with the above mentioned 
feature, the processing means approximately calculates 



14 

the addresses of the projection data to be back 
projected through an interpolation processing based on 
a plurality of representative addresses of detecting 
means in the image data segment region. In other words , 
5 in the back projection processing for every image data 
segment, the addresses of the detecting means on the 
extracted projection data segment are calculated 
through an interpolation by making use of the plurality 
of the limited representative addresses of the 

10 detecting means on the image data segment region. 
Thereby, the address calculation of the detector means 
in the back projection processing can be performed in 
high speed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 Fig.l is a diagram for explaining a fundamental 

constitutional difference between an SDCT and an MDCT, 
Fig. 2 is a diagram for explaining a general idea 
of the back projection image reconstruction method in 
the SDCT, 

20 Fig. 3 is a diagram for explaining a general idea 

of the back projection image reconstruction method in 
the MDCT, 

Fig. 4 is a diagram showing an entire constitution 
of an MDCT representing an embodiment of a tomographic 
25 device according to the present invention. 

Fig. 5 is a diagram for explaining an example of 
dividing an image reconstruction area in the 
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tomographic device according to the present invention. 

Fig. 6 is a diagram for explaining projection data 
segments cut out in correspondence with image data 
segments, which are constituted by dividing the image 
5 reconstruction area. 

Fig. 7 is a diagram showing a cut out processing 
flow of the projection data segments corresponding to 
the image data segments, and 

Fig. 8 is a diagram for explaining a general idea 
10 of an interpolation processing at step S 86 in Fig. 7. 

BEST MODES FOR CARRYING OUT THE INVENTION 

[0019] 

Herein below, an embodiment of the tomographic 
device according to the present invention will be 

15 explained in detail with reference to the drawings 
accompanied. Fig. 4 is a diagram showing an entire 
constitution of an MDCT representing an embodiment of 
a tomographic device according to the present invention. 
A scan method of the MDCT is a rotate-rotate method (the 

20 third generation) and the MDCT, when roughly sectioned, 
is constituted by a scanner 40, an operation unit 50 
and a bed 60 for moving a subject while setting the same 
thereon . 

[0020] 

25 The scanner 40 is constituted such as by a central 

control device 400, an X-ray control device 401, a high 
voltage generating device 402, a high voltage switching 
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unit 403, an X-ray generating device 404, an X-ray 
detector 405, a pre-amplif ier 406, a scanner control 
device 4 07, a scanner drive device 4 08, a collimator 
control device 409, a bed control device 410 and a bed 
5 movement measuring device 411. The operation unit 50 
is constituted by an input-output device 51 including 
such as a display device, an input device and a memory 
device and a computing device 52 including such as an 
image reconstruction computing device and an image 

10 processing device. The input device is constituted by 
such as a mouse and a key board, and is for such as 
measuring bed moving speed information and an image 
reconstruction position and for inputting parameters 
of image reconstruction, the memory device is for 

15 storing these information and the display device is for 
displaying these information and a variety of data such 
as reconstructed images. The image reconstruction 
computing device is for processing projection data 
obtained from a multiple row detectors and the image 

20 processing device is for applying a variety of 
processing to the reconstructed images and for 
displaying the same. 
[0021] 

The central control device 4 00 transmits control 
25 signals necessary for imaging to the X-ray control 
device 4 01, the bed control device 410 and the scanner 
control device 4 07, based on command inputs from the 
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input device in the operation unit 50 with regard to 
such as scanning conditions (such as bed moving speed. 
X-ray tube current. X-ray tube voltage and slice 
position) and reconstruction parameters (such as a 
5 region of interest, reconstructed image size, back 
projection phase width and reconstruction filter 
function) , and begins an imaging operation upon receipt 
of an imaging start signal. When the imaging operation 
begins, a control signal is sent to the high voltage 

10 generating device 402 from the X-ray control device 401, 
a high voltage is applied to the X-ray generating device 
404 via the high voltage switching unit 403, X-rays 
emitted from the X-ray generating device 4 04 are 
irradiated to a subject and the penetrated light thereof 

15 makes incident to the X-ray detector 405. At the same 
time, from the scanner control device 407 a control 
signal is sent to the scanner drive device 4 08 and the 
X-ray generating device 404, the X-ray detector 405 and 
the pre-amplif ier 406 are controlled so as to rotate 

20 circumf erentially around the subject. 
[0022] 

X-rays emitted from the X-ray generating device 
404 are controlled of their irradiation area by a 
collimator 412 controlled by the collimator control 
25 device 4 09, absorbed (attenuated) by respective tissues 
within the subject, penetrate the subject and detected 
by the X-ray detector 405. The X-rays detected by the 
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X-ray detector 4 05 are converted therein into an 
electric current and which is amplified by the 
pre-amplif ier 406 and is input as projection data 
signals to the computing device 52 in the operation unit 
5 50. The projection data signals input to the computing 
device 52 are subjected to an image reconstruction 
processing at the image reconstruction computing device 
in the computing device 52. The reconstructed image is 
stored at the memory device in the input-output device 
10 51 and is displayed as a CT image on the display device 
in the input-output device 51. 
[0023] 

As shown in Fig.l, different from the single row 
detector type CT, in the multiple row detector type CT, 

15 since detector elements are arranged in a plurality of 
rows in the circumferential revolving axis direction, 
as a whole, a detector having a broader width than the 
single row detector type CT is realized. Further, in 
the single row detector CT, the X-ray beams thereof 

20 cross perpendicularly to the circumferential revolving 
axis, on the other hand, in the multiple row detector 
type CT, the X-ray beams show slanting angles (cone 
angle) with respect to the circumferential revolving 
axis as the beams move away from the mid plane (center 

25 row) of the detector rows. 
[0024] 

Fig. 5 is a diagram for explaining an example of 
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dividing an image reconstruction area in the 
tomographic device according to the present invention. 
In the present embodiment^ through a program stored in 
the central control device 400, 3-dimensional image 
5 reconstruction area is divided into image data segments 
ei'-Sp of P pieces = M x N x L. When the matrix of image 
to be reconstructed is 512x512x512, the number of 
division is determined as follows while assuming the 
number of division in x axis direction is M, the number 
10 of division in y axis direction is N and number of 
division in z axis direction is L; 

M=2^ m is an integer of more than 0 
N=2" n is an integer of more than 0 
L=2'^ 1 is an integer of more than 0 
15 [0025] 

By determining the number of division as above, 
the image reconstruction area can be divided in an equal 
integer unit size respectively in x, y and z directions, 
thereby, a processing loop for accessing to respective 
20 pixels within the image data segments ei'-Gp can be 
shared in common through out the image data segments 
ei-^Sp, which reduces processing complexity. 

Fig. 6 is a diagram for explaining projection data 
segments cut out in correspondence with image data 
25 segments, which are constituted by dividing the image 
reconstruction area. In the back projection processing 
according to an embodiment of the present invention. 



20 

at firsts the projection data input in the computing 
device 52 are divided and extracted through a programs 
in the image reconstruction computing device in the 
computing device 52 into projection data segments 
5 of small size, which are necessary for reconstructing 
respective image data segments 61-63 and which are read 
in and stored in the high speed memory. Then, through 
another program stored in the image reconstruction 
computing device in the computing device 52 the back 

10 projection processing is performed based on the 
projection data segments ll'-lS stored in the high speed 
memory. Further., in this embodiment, for the 
simplicity's sake., the projection data regions are 
extracted in a rectangular shape, the regions can be 

15 extracted in a form of a polygon such as rhombus and 
a parallelogram. 
[0026] 

Fig. 7 is a diagram showing a cut out processing 
flow of the projection data segments corresponding to 

20 the image data segments. In calc_address () at step S 
81, respective detector addresses (rwl,chl), 
(rw2,ch2) , (rw3., ch3)., (rw4, ch4) (addresses on detectors 
of X-rays passing through four points at corners.) of 
the projection data segment 71 corresponding to the 

25 four points at the corner 

(p(xl,yl) ,p(x2,y2) ,p(x3,y3) ,p(x4,y4) ) of the 
2-dimensional image data segment 61 are calculated 
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through a program stored in the image reconstruction 
computing device in the computing device 52. Further, 
in case of 3-dimensional image data segment, respective 
detector addresses (rwl, chl.) , . . . (rw8,ch8.) of the 
5 projection data segment corresponding to the eight 
points at the corner 

(p(xl,yl,zl) ,p(x2,y2,z2) , . • .p(x8,y8,z8) ) are 
calculated. For the above calculation, the equations 
(l)-'(6) as mentioned above can be used. 
10 [0027] 

In calc_maxmin4 () at step S 82, respective maximum 
value and minimum value (max_rw., max_ch, min_rw., min_ch) 
among the four detector addresses calculated at step 
S 81 are calculated. These values can be calculated by 
15 simply comparing the values of the detector addresses 

{ rwl , chl ) , ( rw2 , ch2 ) , (.rw3., ch3 ) (.rw4 , ch4.) . 

[0028] 

In calc_cut_size () at step S 83, the sizes 
.(rw_size., ch_size) in the row direction and in the 

20 channel direction of the projection data segment 71 as 
shown in Fig. 6 are calculated based on the detector 
addresses calculated at step S 82. These sizes can be 
calculated by substituting the maximum value and 
minimum value (max_rw, max_ch, min__rw., min_ch) of the 

25 detector addresses calculated at step S 82 for the 
following equations ; 

rw size=max rw-min rw 
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ch_s i z e=max_ch--inin_ch 
[0029] 

In calc_base_address () at step S 84, a reference 
address (rw_base^ ch_base) on the projection data of the 
5 projection data segment 71 as shown in Fig. 6 is 
calculated- This address can be calculated by 
substituting the maximum value and minimum value 
(max_rw,max_ch,min_rw,min_ch) of the detector 
addresses calculated at step S 82 for the following 
10 equations; 

rw_base=min_rw 

ch_base=min_ch 

[0030] 

In cut_data () at step S 85, the projection data 
15 segment 71 is extracted from the projection data based 
on the size {rw_size, ch_size) of the projection data 
segment 71 and the reference address (rw_base, ch_base) . 
In the present embodiment, although the size of the 
projection data segment is calculated at step S 83, a 
20 predetermined fixed size, which is sufficient for 
storing the projection data segment, can be used. 
[0031] 

In calc_interpolation_data () at step S 86, 
detector addresses at reconstruction points in the 
25 image data segment are calculated from a plurality of 
limited but not all detector addresses in the image data 
segment 61. Fig. 8 is a diagram for explaining a general 
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idea of an interpolation processing at step S 86. As 
shown in Fig. 8, through a linear interpolation 
processing with a program stored in the image 
reconstruction computing device in the computing device 
5 52 a detector address of a reconstructing pixel point 
p(x5,y5)is calculated based on the detector addresses 
corresponding to four corner points 

p (xl, yl) ,p (x2, y2) ,p (x3,,y3) ,p (x4, y4) of the divided 
image data segment 61 in a rectangle. Specifically/ when 

10 assuming that the detector addresses corresponding to 
the four corner points 

p(xl,yl) ,,p(x2, y2) ,.p (x3, y3) ,p (x4 , y4 ) as 
( rwl , chl ) , < rW2 , ch2 ) , ( rw3 , ch3 ) , ( rw4 , ch4 ) , the detector 
address of the reconstruction point as (rw5,ch5) is 

15 determined as follows while multiplying respective 
interpolation coefficients to the coordinate values and 
adding the same; 

rw5=coefn'»Twl+coeff2*rw2+coefr3*rw3+coeff4*rw4 
ch5=coefFl*chl+coeff2*ch2+coeff3*ch3+coeflF4*ch4 

20 

[0032] 

Wherein, coef f 1, coef f2, coef f 3,coef f 4 are the 
interpolation coefficients and in the case of Lagrange 
interpolation, these are determined as follows; 

^ ^ coeffl=((x5-x2) * (x5-x3) * (x5-x4)) / ((xl-x2) * (xl-x3) * (xl-x4)) * 

((y5-y2) * (y5-y3) * (y5-y4)) / ((y l-y2) * (yl-y3) * (y l-y4)) 
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coefl2=((x5-xl) * (x5-x3) * (x5-x4)) / ((x2-xl) (x2-x3) ^< (x2-x4)) * 

((y5-y 1) * (y5-y3) * (y5-y4)) / ((y2-y 1) * (y2-y3) * G^2-y4)) 

coeff3=((x5-xl) * (x5-x2) * (x5-x4)) / ((x3-xl) * (x3-x2) * (x3-x4)) * 
5 ((y5-yl) * (y5-y2) * (y5-y4)) / ((y3-yl) * Cv3-y2) * (y3-y4)) 

coeff4=((x5-xl) * (x5-x2) * (x5-x3)) / ((x4-xl) ^< (x4-x2) (x4-x3)) ^< 

((y5-yl) * (y5-y2) * G'5-y3)) / ((y4-y 1) M-y2) * (y4-y3)) 



10 In this embodiment^ although the four corner 

points are used^ the present invention is not limited 
thereto. Further^ a six points interpolation determined 
from up and down and right and left points can be used. 
Still further^ such detector address calculation 

15 through interpolation can be performed only for 
determining row position of which calculation is 
complex. After completing the detector address 
calculation of the reconstruction point in the image 
data segment 61 at step S 8 6, like processing is executed 

20 for the subsequent image data segment 62 and the 
projection data segment 72. At this instance, since two 
points p (x2, y2) ,p (x3, y3) of the image data segment 61 
and the detector addresses (rw2,ch2), <rw3,ch3) of the 
projection data segment 71 can be shared in common, the 

25 processing complexity and memory capacity increase can 
be reduced. Herein above for the sake of simplifying 
explanation, the interpolation processing is primarily 
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explained with reference to 2-dimensional image data, 
like process can be used for 3-dimensional image data. 
[0033] 

Herein above, the embodiment of the present 
5 invention has been explained in detail, however, such 
is simply intended for explanation and exemplification 
and the present invention is not limited thereto. 

Further, in the present embodiment, although an 
example of the tomographic device using X-ray has been 
10 explained, the present invention is not limited thereto 
and the present invention is applicable to a tomographic 
device using such as neutron beams, positron beams, 
gamma beams and light beams. 
[0034] 

15 Further, the scanning method of the present 

invention is not limited to any one of those for the 
first generation, second generation, third generation 
and fourth generation and the present invention is 
applicable to such as a multi tube CT mounting a 

20 plurality of X-ray sources, a cathode scanning CT and 
an electron beam CT. Still further, the present 
invention can be applicable to a variety of detectors 
having different configurations such as a detector 
arranged around a cylindrical surface while disposing 

25 an X-ray source at the center, a flat detector, a 
detector arranged on a spherical surface while 
disposing an X-ray source at the center and a detector 
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arranged around a cylindrical surface while locating 
the circumferential revolving axis at the center. Still 
further, the present invention is not limited to a 
spiral orbit scan, but can be applicable to a circular 
5 orbit scan. Still further, although in the above the 
image reconstruction area was divided in a same number 
in X and y directions., the present invention is not 
limited thereto and the image reconstruction area can 
be divided in different numbers in x and y directions. 

10 Still further, in the above the image reconstruction 
area was divided in a rectangle in x, y and z space, 
the image reconstruction area can be divided in a 
polygon such as a triangle and an octagon, and further, 
the image reconstruction area can be divided in a polar 

15 coordinate. 

[0035] 

Still further, with regard to the division of the 
image data segments to P=M x N x L, the division number 
can be designed to be input through a provision of a 

20 measure which permits such input externally via the 
input-output device 51. Further, the projection data 
input to the computing device 52 can be designed 
displayable through a provision of a measure, which 
permits the display on the display device in the 

25 input-output device 51 together with the divided image 
data segment regions. More specifically, through 
connecting the projection data, projection image from 
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front and/or side of a subject being laid down is 
displayed as scanogram, a 3-dimensional image 
reconstruction area is set on the scanogram and how the 
set 3-dimensional image reconstruction area to be 
divided into the image data segments is designed to be 
displayable such as by drawing division lines in the 
3-dimensional image reconstruction area expressed such 
as a rectangle and a square. Further, it can be designed 
that any image data segments can be selectable from the 
scanogram displayed through a provision of a measure, 
which permits the selection thereof externally via the 
input-output device 51. 



